import pandas as pd

wD = pd.read_excel('测风数据2.0（完整）杨明海&焦志雄 .xlsx',parse_dates=['Date & Time Stamp'])
den = 1.041
#print(wD['Date & Time Stamp'].dt.month)

wD['month']=wD['Date & Time Stamp'].dt.month
monthDay={'1月':31,'2月':28,'3月':31,'4月':30,'5月':31,'6月':30,'7月':31,'8月':31,'9月':30,'10月':31,'11月':30,'12月':31}
monthAvgWV={'1月':0,'2月':0,'3月':0,'4月':0,'5月':0,'6月':0,'7月':0,'8月':0,'9月':0,'10月':0,'11月':0,'12月':0}
monthAvgWPD={'1月':0,'2月':0,'3月':0,'4月':0,'5月':0,'6月':0,'7月':0,'8月':0,'9月':0,'10月':0,'11月':0,'12月':0}

n=0
for month in wD['month']:
    m=str(month)+'月'
    monthAvgWV[m]+=wD['CH1Avg'][n:n+1].values               #其中CH1Avg代表70m高度的数据，其他高度数据只需要对这个数据进行修改即可
    monthAvgWPD[m]+=1/2*den*(wD['CH1Avg'][n:n+1].values)**3
    #print(wD['CH1Avg'][n:n+1])
    n+=1
    #wD['month'][n]=m

for month in monthAvgWV:
    monthAvgWV[month] = format(float(monthAvgWV[month]/monthDay[month]/(24*6)),'.3f')
    monthAvgWPD[month] = format(float(monthAvgWPD[month]/monthDay[month]/(24*6)),'.3f')
    print(month,'平均风速:',monthAvgWV[month])
    print(month,'平均风功率密度:',monthAvgWPD[month])

print(type(monthAvgWV['1月']))
#print(monthAvgWV)

#print(wD['Date & Time Stamp'][2]) 

